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We report results from full QCD calculations with two flavors of dynamical staggered fermions on anisotropic 
lattices. The physical anisotropy as determined from spatial and temporal masses, their corresponding dispersion 
relations, and spatial and temporal Wilson loops is studied as a function of the bare gauge anisotropy and the bare 
velocity of light appearing in the Dirac operator. The anisotropy dependence of staggered fermion flavor symmetry 
breaking is also examined. These results will then be applied to the study of 2-flavor QCD thermodynamics. 



1. INTRODUCTION 

Anisotropic lattices have been used extensively 
for T = calculations 0| . Finite temperature 
calculations using anisotropic lattices exploit the 
natural asymmetry of finite temperature field the- 
ory to reduce lattice spacing errors associated 
with the transfer matrix at less cost than is re- 
quired for the full continuum limit p| . 

With a sufficiently small value for the temporal 
lattice spacing, at, we can vary the temperature 
in small discrete steps by varying the number of 
time slices, N t . In this study we sweep through 
the transition by studying different N t 's (8-64). 
Varying the temperature at fixed temporal and 
spatial lattice spacing separates temperature and 
lattice spacing effects, allowing a study of the 
temperature dependence with all other parame- 
ters fixed. 

The coefficients needed to relate lattice observ- 
ables to the physical energy and pressure are de- 
termined as a by-product of the zero temperature 
studies needed to choose the bare parameters. 
Once determined, these "Karsch" coefficients || 
can be used for all temperatures since they de- 
pend only on the intrinsic lattice parameters and 
not on Nt- This allows a straight-forward deter- 
mination of the temperature dependence of the 
energy and pressure, again at fixed lattice spac- 
ing. With two or more slightly different values 
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for at , a high- resolution sampling of temperatures 
can be investigated. 

As the temporal lattice spacing, a t , approaches 
the continuum limit (a t — > 0), the part of the 
flavor symmetry, which is violated by terms of 
0(af), is expected to be restored. In this study 
we are examining our data for evidence of im- 
provement of the flavor symmetry, when at be- 
comes sufficiently small. 

2. THE ANISOTROPIC STAGGERED 
ACTION 

We are simulating full QCD with two dynami- 
cal flavors of staggered fermions on an anisotropic 
lattice. Our calculations are based on the QCD 
Sp, where the gauge action is: 



action = Sq 
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and the fermion action is: 
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In our simulations we attempt to examine the 
QCD phase transition for volumes (16 3 x4), quark 
masses (m w /m p rs 0.3) and a spatial lattice spac- 
ing (a s k 0.3 fm) similar to those used in Nt =4, 
2-flavor thermodynamic studies on isotropic lat- 
tices. Thus, we adjust the bare bare anisotropy 
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(£o) an d the renormalization of the speed of light 
(y±) to yield the required a s and m^/nip but 
work with a much smaller temporal lattice spac- 
ing choosen so the critical value of N t is approxi- 
mately 16. 



3. SIMULATIONS 



run volume 



traj. (3 



m f 



1 


16 3 x32 


5800 


5.425 


1.5 


0.025 


2 


16 2 x24x32 


5100 


5.425 


1.5 


0.025 


3 


16 2 x24x64 


1300 


5.695 


2.5 


0.025 


4 


16 2 x24x64 


1400 


5.725 


3.44 


0.025 


5 


16 2 x24x64 


3400 


5.6 


3.75 


0.025 


6 


16 2 x24x64 


3200 


5.3 


3.0 


0.008 


7 


16 3 x24 


3500 


5.3 


3.0 


0.008 


8 


16 3 x20 


2800 


5.3 


3.0 


0.008 


9 


16 3 xl6 


4500 


5.3 


3.0 


0.008 


10 


16 3 xl2 


8200 


5.3 


3.0 


0.008 


11 


16 3 x8 


5500 


5.3 


3.0 


0.008 



Table 1 

Parameters of all 
dynamical v t = I 
1.2. 



the calculations. All runs have 
except run 3 which has v t — 



We have used the zero temperature runs 1-6 for 
scale-setting and determination of the anisotropy. 
For all the zero-temperature runs the renormal- 
ized anisotropy £ r is calculated both from the p 
masses in the spatial and temporal directions and 
from matching the static potentials (J]. Figure 1 
shows that both methods provide values for £ r 
which are reasonably close. 

The idea behind run 7 through 11 is that we 
keep the spatial lattice spacing, a s , and all other 
run parameters constant and change only the 
number of lattice points, Nt, in the temporal di- 
rection. 

Our simulations implement the R-algorithmj^] 
with step-size At = 0.005 for all jobs from Ta- 
ble 1. Figure 2 shows that our choice of At is 
consistent with the requirements for small finite 
step-size error on the physical quantities. 



4. THE VELOCITY OF LIGHT 

The velocity of light on a lattice can be defined 
through the meson dispersion relation: 



Et,phys{Ps 



2 1 
' Ps.lat ~2 

at 



We tune v t , so that Ct s (p s ) ~ 1- The velocity 
of light is calculated for the it propagating in the 
temporal direction and having non-zero momen- 
tum for three values of v t : 1.0, 0.8 and 1.2, where 
the first is a dynamical parameter and the last 
two are valence parameters. From Figure 3 we 
see that the choice of v t = 1.0 gives velocity of 
light closest to 1.0. 

5. IMPROVEMENT OF THE FLAVOR 
SYMMETRY 

Table 2 shows the meson masses for run 3 (a s — 
0.231(7) fm, a t = 0.072(4) fm) and run 4 (a, = 
0.243(3) fm, a t = 0.057(2) fm). 

We choose A w = (m T2 — m T )/m p , where 7T2 is 
the second local staggered pion, as a quantitative 
measure of the flavor symmetry breaking in the 
spatial (S) and temporal (T) directions. The data 
shows that in the temporal direction for both runs 
A w is smaller than its value in the spatial direc- 
tion, which means that we are seeing improve- 
ment of the flavor symmetry as a t becomes finer. 
Especially for run 4, the ir and 7r 2 look virtually 
degenerate. 



am 


S, #3 


T, #3 


S, #4 


T, #4 


am-jr 


0.680(3) 


0.214(2) 


0.771(4) 


0.181(2) 


am- K2 


0.759(12) 


0.219(2) 


0.836(21) 


0.181(2) 


am p 


0.902(26) 


0.283(4) 


0.948(13) 


0.222(4) 


A. 


0.09(1) 


0.017(9) 


0.07(2) 


0.000(2) 



Table 2 

Meson masses for run 3 and 4. 



6. THE PHASE TRANSITION 

The thermodynamics runs 7-11 have volume 
16 3 x N t , where N t = 8, 12, 16, 20 and 24. The 
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bare parameters are kept constant for all runs 6— 
11. Figure 4 shows the sweep through the phase 
transition as we gradually change the tempera- 
ture by varying only N t . We fit our data to a 
hyperbolic tangent form and determine T c w 160 
MeV from the inflection point. 

7. CONCLUSIONS 

We have studied the finite temperature QCD 
phase transition using staggered fermions on an 
anisotropic lattice with anisotropy of « 4. This 
allows us to explore the temperature dependence 
of the transition with all other parameters fixed. 
The results are roughly consistent with earlier, 
isotropic studies, showing a value of the critical 
temperature of approximately 160 MeV. While 
this approach naturally reduces finite lattice spac- 
ing errors associated with at, we plan to include 
improvements to the spatial parts of the staggered 
fermion action so that the 0{a 2 s ) errors are re- 
duced as well. 
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Figure 1. Scatter plots of the zero temperature 
runs 1-6. 
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Figure 2. Step-size test. Volume 8 3 x 32, (3 = 
5.35, v t = 1.0, <; = 3.5, m f = 0.006. 




Figure 3. Tuning of the velocity of light using a 
16 2 x 24 x 64 volume, (3 = 5.3, Co = 3.0, m f = 
0.008. 




Figure 4. The temperature dependence of (ipip) 
in the region of T c . 



